Aging Exponents in Self-Organized Criticality 
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In a recent Letter [Phys. Rev. Lett. 79, 889 (1997)] we have demonstrated that the avalanches 
in the Bak-Sneppen model display aging behavior similar to glassy systems. Numerical results 
for temporal correlations show a broad distribution with two distinct regimes separated by a time 
scale which is related to the age of the avalanche. This dynamical breaking of time-translational 
invariance results in a new critical exponent, r. Here we present results for r from extensive numerical 
simulations of self-organized critical models in d = 1 and 2. We find ra=i = 0.45 ± 0.05 and 
rd-2 ~ 0.23 ± 0.05 for the Bak-Sneppen model, and our results suggest r — 1/4 for the analytically 
tractable multi-trade model in both dimensions. 



PACS numbers: 64.60.Lx, 05.40. +j, 05.70.Ln 



I. INTRODUCTION 



Self-organized criticality (SOC) describes a general 
property of slowly-driven dissipative systems with many 
degrees of freedom to evolve intermittently in terms of 
bursts spanning all scales up to the system size. Many 
natural avalanche-like phenomena have been represented 
using this concept, including earthquakes [|]-f3, extinc- 
tion events in biological evolution and landscape 
formation [||,||. Recently, SOC has been observed in con- 
trolled laboratory experiments on rice piles pToj ] . Theo- 
retical models of rice piles Q are related to a variety 
of different physical systems by universality Jl2| . One 
crucial ingredient for a system to exhibit SOC is the ex- 
istence of thresholds that allow it to record the stress 
exerted by the driving force over long periods of time. 
The emergence of long-term memory has been demon- 
strated analytically Q for a multi-trait evolution model 
fl4|| , a variant of the Bak-Sneppen model 0. 

In a recent Letter Jl5| we have shown that the self- 
organized critical state in the Bak-Sneppen model ex- 
hibits aging behavior that is reminiscent of glassy sys- 
tems [ ^6[ . Our results indicate that intrinsic two-time 
autocorrelation functions P(t w ;t), describing the return 
of activity to a site at time t a = t + t w which was ac- 
tive most recently at time t a — t w for an avalanche that 
started at t a = 0, decay as power laws with two distinct 
regimes according to 



constant (x <C 1), 

x- r 



(1) 



The early time regime is that of the familiar station- 
ary dynamics. The late time regime has a new critical 



coefficient r characterizing the nonstationary relaxation 
behavior of the SOC systems. The "waiting time" t w 
separating the early and late time regimes is a measure 
of the age of the avalanche. We have argued that this 
aging behavior arises from the hierarchical structure of 
the avalanches. 

Generally, the origin of aging behavior does not have 
to be profound. For instance, it arises in a simple random 
walk model near a wall, where a symmetry (translational 
invariance) is explicitly broken. But while the random 
walk near a wall and the Bak-Sneppen model (especially 
its random neighbor variant fl7|| ) are similar in many 
ways, it appears that those similarities do not explain 
the observed aging behavior. 

In Sec. II, we use the random walk near a wall to illus- 
trate the quantities that will be measured for the SOC 
models. In Sec. Ill, we consider the important technical 
issue of extracting the intrinsic aging behavior for a pro- 
cess that does not conserve its norm (such as a random 
walk near an absorbing wall or a SOC avalanche, having 
a finite stopping probability). In Sec. IV, we present 
detailed numerical results for the aging behavior of SOC 
models in one and two dimensions. In particular, we have 
simulated the Bak-Sneppen model and the multi-trade 
model. We have also simulated sandpile models exhibit- 
ing SOC which show a quite different behavior and will 
be discussed elsewhere In Section V we discuss our 
results and show that a simple random walk description 
is not sufficient to explain the found aging behavior. 



II. AGING RANDOM WALKS 



In this section we will discuss a random walk model to 
provide a simple intuitive picture of aging. It shows how 
the system memorizes the (here, explicit) breaking of a 



1 



symmetry. Furthermore, the random walk illustrates the 
meaning of the correlation functions used to describe the 
aging behavior and serves to discuss some of the technical 
issues in measuring those functions. 

We want to consider a random walker on a d = 1 lattice 
who can jump at most one step on each update either on 
the infinite lattice or a semi-infinite lattice with an ab- 
sorbing wall at the origin. A random walk is completely 
described by its propagator, the conditional probability 
G(n, t\no, to) for a walker to reach a site n at time t, given 
that it was at site no at some previous time to < t. To 
determine its aging behavior we want to compute a sim- 
ple two-time correlation function (see e. g. Ref. jl9| for 
a similar definition) 

P(t w ;t) =J2 G(n,t + t w \n, t w ) G(n, t w \n , 0) (2) 



for a walker to return to a site at time t a = t + t w , given 
that it was at the same site at the "waiting time" t w af- 
ter the start of the walk. Thus, to determine the return 
probability to a site, we take the time t w from the start of 
the walk to a previous passage at that site into account. 
If the two-time correlation function P explicitly depends 
on t w , the system is said to "age" because the walk would 
retain a memory of the time since its inception in form 
of a return probability that evolves over time. 



A. Unconstrained Random Walks 



Clearly, for an unconstraint walk, G is invariant in 
space and time, and it is G(n, t + t w \n, t w ) = G(0, i|0, 0). 
Since the norm of the walk is preserved at all times, 
J2 n G(n,t w \n ,0) = 1, Eq. (|) gives P(t w ;t) = P(t) = 
G(0, t\0, 0), independent of t w . Thus the unconstrained 
walk has no memory of its past and does not age. 



B. Random Walks near an Absorbing Wall 



In the presense of an absorbing wall at the origin, spa- 
tial invariance is explicitly broken while time invariance 
for G still holds. Eq. (0) merely simplifies to 



P(t w ;t) =Y^G(n,t\n,0) G(n,t w \n o ,0) 



(3) 



n>0 



which remains dependent t w . Hence the breaking of spa- 
tial invariance in G leads to a breaking of time invari- 
ance in P. Such a memory effect arises in the follow- 
ing manner: While in the unconstrained case the mean 



distance < n > of the walker from its origin vanishes, 
the walker starting near a wall departs from it such that 
< n >^ t 1 / 2 , which follows from the propagator given 
by 



G(n,t\n ,0) 



1 



u — 



(" + "0)" 



(4) 



for sufficiently large t and n. Since the distribution for all 
walks is sharply peaked near its mean, most walkers oc- 

1/2 

cupy sites n ~ tw away from the wall after the waiting 
time t w . Given that most walkers occupy such sites n, re- 
turns to that site during subsequent times t <^ t w follow 
the statistics of an unconstrained walk: t is not yet large 
enough to return to the wall. Only after times t ~ t w do 
a sizable portion of the walkers experience the effect of 
the distant wall again, which leads to a change in their 
return statistics for all times t^$>t w . The change-over in 
its return behavior at later times t + t w thus provides the 
walker with a memory of the earlier period from the start 
up to time t w when the walk was drifting away from the 
wall. (Considering the walker as the center of a growing 
domain and its distance to the wall as a measure of the 
linear size of that domain relates this random walk model 
nicely to the domain-growth picture of aging in glasses 
proposed in Ref. [p0| .) 

For this random walk model the crossover in its return 
behavior can be easily derived explicitly: Inserting the 
appropriate forms of G in Eq. (jj) into Eq. (^) and choos- 
ing no arbitrarily close to the wall, we obtain asymptot- 
ically 



P(t v 




i 



(5) 



Thus, time-translational invariance is broken and the 
walk appears to age because P becomes a function of 
the scaling variable x = t/t w signaling the predicted 
crossover in the return behavior that scales linearly with 
the waiting time t w . Here the scaling function f(t/t w ) 
behaves asymptotically as f{x <C 1) ~ 1 and f(x ^> 
1) ~ 4/x, but to obtain the "intrinsic" aging behavior of 
this process, we have to consider the effect that the norm 
is not preserved because at each time step walkers may 
disappear at the absorbing wall. 



III. INTRINSIC VS. MEASURED AGING 
BEHAVIOR 



Consider the unconstrained random walk in Sec. [I A, 
but with a finite probability z to disappear at each 
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time step. Then, the propagator G z of this walk is 
given by G z (n, t | no, to) = (1 — z) l ~ to G{n, t\n , to), where 
G is the propagator of the norm-preserving walk in 
Sec. OA. Thus, with J^n G{n, t\no, t ) = 1, we get 
Yl n Gz{n, t\n , t ) = (1 — i. e. the norm of this 

process is not preserved. According to our definition of 
the two-time auto-correlation function P in Eq. (Q) 

P mcas (t w ;t) = ^2G„(n,t + t w \n,t V) ) G z (n,t w \n ,0) 

n 

= G(0,t|0,0)(l-z) t+t ™, (6) 

we would have to conclude that this process ages, since 
P depends on t w . And, indeed, in a numerical simula- 
tion of the process we would measure p mcas j because we 
would average over all processes up to a temporal cut-off 
t co , including those that disappear at times t < t co . 

But, clearly, the "aging" in this simple process is an 
artifact due to the diminished norm (1 - z ) t+t ™ at time 
t + t w . Thus, proper normalization is required to extract 
the "intrinsic" aging behavior (due to the infinite walk) 
from the "measured" aging behavior ^jj . 

Other process, such as random walks near an absorbing 
wall or avalanches in the Bak-Sneppen mechanism below, 
also may disappear before reaching the cut-off, and we 
have to consider the effect on the statistics of the mea- 
sured results. Here, too, one is interested in the intrinsic 
properties of the surviving process, i. e. those of the 
infinite random walk or avalanche, while in simulations 
one usually averages over all runs of a process, whether 
they survived or not. In many cases, the results for the 
asymptotic scaling behavior of some intrinsic property 
are no different than the measured ones because the con- 
tribution from dying runs remains insignificant. But for 
the correlation functions considered in this paper which 
depend on two independent time variables, we do need 
to consider the effect of the probability Pt{0) of a pro- 
cess to disappear at time 9 to relate the intrinsic and the 
measure results. 

To obtain the intrinsic properties of the infinite pro- 
cess, we have to properly normalize the correlation func- 
tion. To that end, we consider the two-time correlation 
function P(t w ;t\9) for a run that disappears exactly at 
time t a = 9, and its generic relation to the intrinsic two- 
time correlation function P lntI (t w ;t): 

fO (9<t w +t), 

P(t w ;t9) = \ (7) 
{P int *(t w ;t) (9>t w +t). 

These quantities are related to the measured two-time 
correlation function P = p mcas given in Eq. (||) : Assum- 
ing a power-law probability Pt{9) ~ 9~ T , r > 1, for the 
run to disappear at time t a = 9, we have 



P m ™ s (t w ;t)= / d9P(t w ;t\9) P t {9) 
Jo 



l-T ,1-T 



p int \t w -t) (t w +ty- T -t 



(8) 



Assuming that we only consider data sufficiently far from 
the cut-off, i. e. (t w + t) ~ r <C t\~ T , we obtain 



P intI (t w ;t)~P mcas (t w ;t) (t w +ty 



(9) 



For the particular form of the intrinsic two-time corre- 
lation function considered in Eq. ([!]), the correct scaling 
function for the aging behavior of the process is given by 



/ intr (x) ~ r cas (x) (i + x y 



(10) 



Since we are interested in the intrinsic behavior f lntl '(x ^ 
1) ~ x~ r , we obtain from our numerical data 



/ mcas (a;) ~ x- [r+T - 1] . 



(11) 



Of course, for the random walk near the absorbing 
wall it is r = 3/2 from the familiar first-passage time 
p2^| . Thus, even after correcting for the effect of dis- 
appearing walkers, the intrinsic process still shows aging 
behavior: The measurable aging effect derived in Eq. (|^), 
f(x) = f mcas (x) ~ a;™ 1 , leads to an intrinsic aging be- 
havior of / lntr (x ^> 1) ~ x~ r with r = 1/2 according to 
Eq. ([ll]) . Interestingly, the infinite walk may never touch 
the absorbing potential on the wall but clearly feels its 
effect. 



IV. AGING IN SELF-ORGANIZED CRITICAL 
MODELS 



The Bak-Sneppen model has been studied intensely 
and with great numerical accuracy in recent years. We 
refer to Ref. |23| for a review of its many features and 
simply utilize those facts here. The model consists of 
random numbers Xi between and 1, each occupying a 
site i on a d-dimensional lattice. At each update step, the 
smallest random number \ m in(t) is located. That site 



as well as its 2d nearest neighbors each get new random 
numbers drawn independently from a flat distribution be- 
tween zero and one. The system evolves to a SOC state 
where almost all numbers have values above A c , with A c 
avalanches formed by the remaining numbers below. 



The multi-trade model |FJ,|13| is a variant of the Bak- 
Sneppen model that provides a series of exact results for 
the spatio-temporal correlations in the avalanche process. 



3 



Especially, an equation of motion can be derived and 
solved to obtain a propagator for the spread of avalanche 
activity, and to obtain a complete set of scaling exponents 
that verify scaling relations previously proposed for the 
Bak-Sneppen model. In this model each lattice site is 
occupied by M independent random numbers. On each 
update again the smallest number in the whole system is 
located and updated, and one randomly chosen number 
(out of M) from each of the neighboring sites is updated 
as well. While the mechanism proceeds in the same way 
for all M as for M = 1 (the Bak-Sneppen model), it can 
be treated analytically for M = oo. 

A quantity similar to the two-time autocorrelation 
function defined in Eq. (||) can be measured for 
avalanches in the self-organized critical state of both 
models. We focus on a simple quantity, Pfi rs t(i), mea- 
suring the first returns of the activity to a given site. 
A power law distribution for P& rs t(t) has been measured 
numerically for a variety of different SOC models |2^] 
by recording all first returns and has been derived ex- 
actly for the multi-trade model [Q. Here we determine 
the intrinsic probability Pfi rs t(itu; t) to return after t time 
steps to a site that was visited most recently at time 
t w from the beginning of the avalanche. Thus, to ob- 
tain the first-return probability, we take the age of the 
avalanche, t w , into account. While in the stationary state 
of SOC models the first return distribution is generally 
•Pfirst(*) ~ i~ rfirBt , (t — > oo), we find that P first (t w ;t) for 
both models considered here scales according to Eq. (|l|) 
where the exponent Tfi rst can be related to other critical 
exponents via scaling relations for SOC f23| . The ori- 
gin of the intrinsic aging exponent r appears to us to be 
non-trivial, signaling the breaking of time-translational 
invariance in the avalanche dynamics. Unlike in the ran- 
dom walk near a wall, no symmetry is explicitly broken. 
The Bak-Sneppen mechanism in both models evolves on 
an isotropic lattice with update rules that do not change 
with time. The question then arises whether the expo- 
nent r in Eq. (I can be related to the known universal 
coefficients of the stationary SOC process, or whether it 
describes new physics in avalanche dynamics of the Bak- 
Sneppen model. 



A. Numerical Procedure 



In our simulations for both models we have used the 
equivalent branching process p4] to eliminate any finite- 
size effects. Initially, at time t a = 0, the smallest thresh- 
old value is set equal to A c to start a A c avalanche. In 
every update t a — > t a + 1, only the signal A m i n (i a ) and 
its 2d nearest-neighbor sites receive new threshold val- 
ues. At any time, we store only those threshold values 



Xi < A c that are part of the avalanche because only those 
numbers can contribute to the signal, i. e. can ever be- 
come smallest number. In addition, we keep a dynamic 
list of every site that has ever held the signal at some 
time to determine the first-return probabilities. (Since 
the Bak-Sneppen mechanism for d < 4 is a fractal re- 
newal process, activity always returns to a site unless 
the avalanche dies.) Avalanches die when there are no 
Xi < A c or are stopped at a cut-off t a = t co , and a new 
(independent) avalanche is initiated with t a = 0. 

At each update t a we determine the previous time t w 
when the signal was on the same site most recently (if 
ever). Then its first-return time is give by t = t a — t w , 
and we bin histograms labeled by i = [| log 2 t w ~\ and 
j = [log 2 t] . The data is binned logarithmically so that 
in each bin a comparable number of events is averaged 
over: for each increment of j, the width of the bins for t 
increases by a factor of 2, while for each increment of i 
the t w bins increase by a factor of 8. We then normalize 
for each value of i separately to obtain the measured first- 
return probabilities P^.^ s (t w ;t). This data is plotted in 
Figs. HI and ^ 

In each Figure, each graph refers to a different value of 
i, increasing by a factor of 8 each time from left to right. 
Each graph possesses two distinct power law regimes, sep- 
arated by a crossover. To determine the form of the scal- 
ing function f(x) for these graphs according to Eqs. ([j]), 
we note that the crossover appears to scale linearly with 
t w in all cases. Thus, we plot 

/""W ~ t T ^P£X{t w ;t), x=± (12) 

using the appropriate values of Tfi rst . In each case, the 
data collapses reasonably well onto a single curve, f(x), 
which is constant for small argument, and appears to fall 
like a power law, see Figs. |H and | The exponent 
of the power law is given by r + t — 1 according to the 
relation between measured and intrinsic data discussed 
in Sec. Ill [see Eq. (fill)]. The values of r are given in 
Ref. (23) for the Bak-Sneppen model in d = 1 and 2, and 
r = 3/2 in any dimension for the multi-trade model. 



B. Results for the Bak-Sneppen Model in d — 1 

We have simulated the Bak-Sneppen branching pro- 
cess in d = 1 with A c = 0.66702, summing over a se- 
quence of all avalanches up to a cut-off at t co = 2 27 . That 
data consists of a total of about 10 11 updates. (The re- 
sults reported here are consistent with but substantially 
better than those reported previously in Ref. M ] where 
data from avalanches that did not reach the cut-off was 
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discarded to avoid confusion about the relation between 
measured and intrinsic properties discussed in Sec. III.) 




FIG. 1. Plot of Pg£?"(t w ;t) as a function of t for the 
Bak-Sneppen model in d = 1. Each graph is offset by a factor 
to avoid overlaps. 




itr 2 kt 1 io° 10 1 io £ 10 3 

x=t/t„ 



FIG. 2. Scalingplot for / mcas (x) as a function of x = t/t w 
according to Eq. ( Jl 2] ) for the normalized distribution of the 
Id Bak-Sneppen model in Fig. [j]. 

The distributions for P^.^ s (t w ; t) as a function of t, 
normalized for each value of i, 8 l_1 < t w < 8 Z , are plot- 
ted in Fig. [Ij for i — 1, . . . , 8. Each graph shows two 
scaling regimes separated by a crossover that appears to 
scale linearly with the associated value of t w . The initial 



regime scales with the familiar exponent rg rs t = 1.58 p3] 
(as indicated by the dashed line to the right). Cut-off 
effects become apparent at about t w 10 7 = 10%t co . 

In Fig. H we have combined the data into a scaling 
plot for / mcas (a:) according to Eq. jl2] ) as a function of 
the scaling variable x — t/t w . For x < 1 that data in- 
deed collapses onto a constant, while we observe a col- 
lapse onto a power law over three orders of magnitude for 
x > 1. Deviations from this behavior is generally due to 
short-time, transient behavior for x < 1, and due to sta- 
tistical noise deep in the tail of the distribution for x > 1. 
We have bracketed the power-law tail by two dashed lines 
~ £-0.55 anc j ^ X -0 - 5 . Thus, we estimate that the expo- 
nent of the tail is given by r + r - 1 = 0.52 ± 0.04. With 
r = 1.07 ± 0.01 [||, we finally get r = 0.45 ± 0.05. 



C. Results for the Bak-Sneppen Model in d — 2 



In this case we have simulated the with A c = 0.328855 
p3[ , summing over a sequence of all avalanches up to a 
cut-off at t co = 2 25 (longer avalanche are less common 
here than in the d = 1 case). We have run about 3 x 10 10 
updates for this model. 

The distributions for P^^ s (t w ;t) as a function of £, 
again normalized for each value of i, 8 4 " 1 < t w < 8\ 
are plotted in Fig. S for i = 1, . . . , 7. As in the d = 1 
case, each graph shows two scaling regimes separated 
by a crossover that appears to scale linearly with the 
associated value of t w . The initial regime supposed to 
scale with the exponent Tfi rst = 1.28 determined from a 
more extensive simulation in Ref. ]23)] . That behavior is 
given by the dashed line to the right. But each graph 
approaches that asymptotic behavior in its initial scaling 
regime only very slowly, indicating strong corrections to 
scaling in this case. Even the combined data (the dashed 
line with square marks), not unlike the corresponding 
plot in Ref. ]23] | , approaches asymptotia very reluctantly. 
Furthermore, cut-off effects become apparent at about 
t w 10 6 = 10%t co . 

In Fig. [| we have combined the data into a scaling 
plot for / mcas (a;) according to Eq. ( |T^ ) as a function of 
the scaling variable x — t/t w . The collapse onto a con- 
stant for the data at x < 1 only proceeds slowly due 
to the aforementioned corrections to scaling in the first 
return probability. On the other hand, we observe a col- 
lapse onto a power law over more than three orders of 
magnitude for x > 1 in this case: Since Tfj rst ~ 1.28 is 
smaller than in the case d = 1, many more events occur 
in the tail of the distribution. We have bracketed the 
power-law tail by two dashed lines ~ x~~ - 5 and ~ x~ 0A5 . 
Thus, we estimate that the exponent of the tail is given 
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by r + t - 1 = 0.47 ± 0.04. With r = 1.245 ± 0.010 [g| 
we finally get about r = 0.23 ± 0.05. 



exponent for the avalanche dimension (D = 4), can be 
derived exactly for this model 14 But many other 
properties are as of yet elusive or can only be inferred 
from scaling relations; for instance, it is 



Tfirst = 2 - 



(13) 
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FIG. 3. Plot of PfosT(i»;i) as a function of t for the 
Bak-Sneppen model in d = 2. Each graph is offset by a factor 
to avoid overlaps. 




FIG. 4. Scaling plot for / mcas (x) as a function of x = t/t w 
according to Eq. ( |l2j ) for the normalized distribution of the 
2d Bak-Sneppen model in Fig. || 



D. Results for the Multi- Trade Model in d = 1 



A number of important properties, like the distribu- 
tions of the width and duration of an avalanche and the 



In particular, no expression for P& T at(tw',t) has been ob- 
tained so far to study the aging behavior in this model 
more explicitly. Thus, as a first step to obtain an explicit 
expression, we have simulated this model also in d = 1 
and 2. The results suggest that r = 1/4 in both cases, 
independent of dimension. 

We have simulated the with A c = 1/2 since we chose 
to update only one number on each neighbor and always 
replace the minimum itself with 1 pfl . Summing over a 
sequence of all avalanches up to a cut-off at t C( 
have run about 10 12 updates for this model. 



)27 



we 



The distributions for P^^ s (t w ; t) as a function of t, 



< 



are 



again normalized for each value of i, 8 < t u 
plotted in Fig. [| for i = 1, . . . , 8. Each graph shows two 
scaling regimes separated by a crossover that appears to 
scale linearly with the associated value of t w . Accord- 
ing to Eq. ([l3|), the initial regime supposed to scale with 
the exponent Tfi rs t = 7/4. That behavior is given by the 
dashed line to the right. But as in the case of the 2d Bak- 
sneppen model we observe strong corrections to scaling. 
Even worse, cut-off effects already become apparent at 
about t w 10 6 = l%t co . Generally, despite of spend- 
ing much more time on the simulation, the data has the 
poorest quality in this case also because the large value of 
rfi rs t suppresses the occurrence of events in the long-time 
tail. 

Nontheless, in Fig. ||, we have combined the data into 
a scaling plot for / mcas (a;) according to Eq. ( [l2"|) as a 
function of the scaling variable x = t/t w . As expected, 
the corrections to scaling prevent a satisfactory collapse 
onto a constant for the data at x -C 1. Apparently, even 
the crossover regime is beset by transient behavior which 
makes it difficult to localize the transition to the power- 
law regime. At best, we can discern scaling over two 
orders of magnitude in the tail before the data gets too 
noisy. On the other hand, it is fair to assume from the an- 
alytical results that any exponent in this model should 
be a multiple of 1/4. Considering that the tail clearly 
scales very close to the dashed line with x~ 0J5 below 
and definitely not like the dashed line with x~ 05 above, 
we believe that r + r-1 = 3/4. Thus, with r = 3/2 @, 
we conjecture r = 1/4. 



G 



FIG. 5. Plot of P£tt B (tw;t) as a function of t for the 
multi-trade model in d = 1. Each graph is offset by a factor 
to avoid overlaps. 




FIG. 6. Scalingplot for / mcas (x) as a function of x = t/t w 
according to Eq. (|12| ) f° r the normalized distribution of the 
Id multi-trade model in Fig. 0. 



E. Results for the Multi- Trade Model in d = 2 



We have simulated the with A c = 1/4 since we chose 
to update only one number on each neighbor and always 
replace the minimum itself with 1. Summing over a se- 
quence of all avalanches up to a cut-off at t co = 2 27 , we 
have run about 10 11 updates for this model. 



FIG. 7. Plot of P§£?"(t w ;t) as a function of t for the 
multi-trade model in d — 2. Each graph is offset by a factor 
to avoid overlaps. 




10 1(T 4 1CT 3 1(T 2 10"' 10° 10 1 10 2 1 3 1 4 1 5 
x=t/t„ 

FIG. 8. Scalingplot for / moas (a;) as a function of x = t/t w 
according to Eq. (|l2j) for the normalized distribution of the 
2d multi-trade model in Fig. [jj. 

The distributions for P^.^ s (t w ;t) as a function of t, 
again normalized for each value of i, 8 l_1 < t w < 8 l , are 
plotted in Fig. ^ for i = 1, . . . , 9. Each graph shows two 
scaling regimes separated by a crossover that appears to 
scale linearly with the associated value of t w . Accord- 
ing to Eq. (|r3|), the initial regime supposed to scale with 
the exponent Tg rs t = 3/2. That behavior is given by the 
dashed line to the right. We observe again strong cor- 
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rections to scaling and cut-off effects already at about 
t w 10 6 = l%t co . 

Combining the data into a scaling plot again, see Fig. ||, 
the corrections to scaling again prevent a satisfactory col- 
lapse onto a constant for the data at x -C 1. But we 
observe an excellent collapse of the data in the tail over 
roughly three to four orders of magnitude, clearly closer 
to the dashed with x -0 - 75 below than the dashed line 
with a: -0,5 above. Thus we have r + t — 1 = 3/4 and 
r = 3/2 as well in this case, leading us once more to 
conjecture r = 1/4. 



V. CONCLUSIONS 



In conclusion, our numerical simulations affirm the ex- 
istence of a new and as-of-yet unexplained power law 
regime in the late time behavior of the Bak-Sneppen 
model that was discover in Ref. The results under- 
score the prevalence of memory effects in the SOC state 
of this model fl§|2§ and possibly other SOC models [@. 

We have been able to rule out a simple relation of the 
new exponent r to the known exponents by considering 
some of the more obvious scaling arguments (which may 
not exist at all) fl5|| . For instance, we can consider an 
avalanche as a random walk near a wall in an abstract 
random number space fl7| . The addition or elimination 
of a random number below the threshold A c corresponds 
to taking a step away or towards an absorbing wall in 
a random walk. The avalanche ends when no random 
numbers are left below A c , i. e. the walker has been 
absorbed at the wall. The growth of random numbers 
below A c scales like < n >~ t ds for the (intrinsic) infi- 
nite avalanche, with d s = 0.11 and 0.25 for the Id and 
2d Bak-Sneppen model [^3) and d s — 1/2 for the multi- 
trade model in any dimension. If the aging behavior in 
these models was due to same effect as in our simple ran- 
dom walk model in Sec. II, we would expect that r = d s 
which is inconsistent with the numerical results. It ap- 
pears that the origin of the aging behavior is due to more 
subtle features of the process and that time-translational 
invariance might be broken dynamically, as we have ar- 
gued in Ref. |lj| . The fact that the analytically tractable 
multi-trade model also shows nontrivial aging behavior 
gives us hope that we will be able eventually to under- 
stand its origin. 

I am very grateful to Maya Paczuski for discussing the 
results in this paper with me. 
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